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The motion of a domain wall in a two dimensional medium is studied taking into account the 
internal elastic degrees of freedom of the wall and geometrical pinning produced both by holes and 
sample boundaries. This study is used to analyze the geometrical conditions needed for optimizing 
crossed ratchet effects in periodic rectangular arrays of asymmetric holes, recently observed experi- 
mentally in patterned ferromagnetic films. Geometrical calculations and numerical simulations have 
been used to obtain the anisotropic critical fields for depinning flat and kinked walls in rectangular 
arrays of triangles. The aim is to show with a generic elastic model for interfaces how to build a 
rectifier able to display crossed ratchet effects or effective potential landscapes for controlling the 
motion of interfaces or invasion fronts. 

PACS numbers: 64.60.Ht,87.85.Qr,75.60.Ch,47.61.-k 



I. INTRODUCTION 



The dynamics of elastic interfaces is responsible for a 
wide variety of physical phenomena in very different ex- 
perimental systems. Prominent examples are the prop- 
agation of reaction fronts or surface growth in mate- 
rial science — , cell motility and membrane dynamics in 
biology^, domain walls in ferromagnetic"^"^ or ferroelec- 
tric films^"— , fluid invasion in porous mediai^, contact 
lines of liquids menisci^, and crack propag 

ationl4'i5_ 

all these cases, the presence of heterogeneities, which lo- 
cally promote wandering, compete with the elasticity of 
the interface, giving rise to complex collective pinning ef- 
fects. Understanding these effects is a challenging prob- 
lem relevant both for the basic and applied viewpoint. 

A particularly interesting case of interface pinning is 
the "geometrical pinning" induced by the presence of ar- 
tificially introduced holes or antidotsi^""— , or by a spatial 
modulation of the sample boundary conditions in nar- 
row samplesiii^r^. These kind of boundaries can pin 
the interface by locally reducing its extension, thus sav- 
ing surface tension energy. For extended domain walls 
this kind of pinning has been recently realized experimen- 
tally and showed to be able to modify the magnetization 
dynamics^ d^i^^" — and to produce, in particular, inter- 



esting ratchet transport of magnetic domain walls^^. 
Being mostly geometrical (i.e. determined mostly by the 
shape and distribution of holes, or by the geometry of the 
boundaries and not much on the specific microscopic pin- 
ning interaction) this kind of pinning has the advantage 
over other artificial pinning mechanisms that it can be 
more easily tailored at a wide range of scales to control 
the wall motion in various specific ways. 

We have recently analyzed, specifically, the pinning ef- 
fect of asymmetric holes on the propagation of domain 
walls in magnetic films, finding that, under certain ge- 
ometries and oscillating external magnetic fields, the mo- 
tion of fiat and kinked walls is rectified in opposite direc- 
tions'^: the asymmetry between forward/backward flat 
wall propagation results in a direct ratchet effect whereas 
the asymmetry between upward/downward kink propa- 
gation along a wall induces an inverted ratchet effect. 
This striking sensitivity yields new strategies to control 
the motion of the wall. The crossed rectification reported 
in Refi^ relies on the difference between the critical fields 
to depin the wall in each direction, an it is also present 
in a generic model for elastic interfaces: the cj)'^ modeP— . 
In this paper we calculate the depinning field of a generic 
(f)'^ interface in the presence of an array of triangular an- 
tidots both by geometrical considerations and numerical 
simulations. We use this simple model because it is the 
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minimum model that captures the essential physics be- 
hind the crossed ratchet effects reported in^. In addition, 
our method can be widely used to design interface rec- 
tifiers of elastic interfaces by using holes or boundary 
conditions in an arbitrary geometry. 

Our starting point is the overdamped 0^ model in the 
plane, i.e., a scalar field (t>{x,y;t) obeying the following 
evolution equation: 



,yc)t0 = cV2(/) + eo(0-0^) + i/ 



(1) 



where c is the elastic stiffness of the order parameter, eg is 
proportional to the local barrier separating two minima 
of the local free energy, H is an external field biasing 
one of the two minima, and the friction coefficient rj sets 
the microscopic time scale. The evolution equation ^ 
derives from the energy functional: 



dxdy 



U{cj>{x,y)) - Hcj>{x,y) + -\V<l>{x,v)\^ 



(2) 

with U{(l)) = eQ{(j)'^ - 1)2/4. 

For H — 0^ and with the appropriate boundary condi- 
tion, say, (j){—L,y;t) = 1; (j){L,y;t) ~ —1, the stationary 
solution of Eq. ([T]) is given by a domain on the left side of 
the plane with positive and approximately homogenous 
field and a domain on the right side with negative field, 
both separated by an interface of width proportional to 
y^c/eQ. When the field is switched on to a positive (neg- 
ative) value, the interface is pushed to the right (left) 
to minimize the total energy. However, the interface has 
also an elastic energy proportional to its length. There- 
fore, if the geometry where the field is defined is such 
that the length of the interface increases when moving to 
the left or right, then the interface will be pinned until 
the field reaches a critical value. Our goal is to provide 
an estimation of such depinning field in a general geom- 
etry and to analyze the geometrical conditions in which 
rectification effects appear in the elastic interface propa- 
gation. 

The organization of the paper is as follows. In Sec. |lTl 
we develop a general theoretical framework to address 
the problem of rectification of domain walls. We reduce 
the field (f> to an elastic wall and derive an analytical ex- 
pression for the local depinning field in two dimensional 
stripes with arbitrarily shaped borders. In Sec. lIIII we ap- 
ply the previous results to build a 2D array of triangular 
holes that displays ratchet effects and give specific pre- 
dictions for the appearance of normal and crossed ratchet 
effects. In Sec. llVlwe summarize our results. 



II. DEPINNING FIELDS IN ARBITRARY 
GEOMETRY: GENERAL THEORY 

A. Prom field equations to elastic walls 

Our aim is to calculate the depinning field of certain 
interfaces in an arbitrary geometry. To simplify the task. 



we first need to reduce the whole field equation ([T]) to 
a parametric description of the interface, in the spirit of 
the collective coordinate approach widely used in one- 
dimensional models^^. 

In Appendix Elwe construct a solution of Eq. ([ij where 
two domains of positive and negative magnetization are 
separated by a wall defined by the line {x{s), y{s)). The 
solution reads: 



4>{x, y) = tanh 



(3) 



where g{x,y) is the distance of point {x^y) to the back- 
bone of the wall (a;(s), y(s)) and 



(4) 



can be considered as its width. Eq. ([3]) is well known as 
the field corresponding to a single wall in the one dimen- 
sional (fA model, and, as we show in Appendix [Aj can be 
extended to two dimensions if the curvature radius of the 
wall {x{s)^y{s)) is much larger than its width w. 

The energy of this solution, for small width w, can be 
approximated by (see Appendix |A| 



E^ad- 2HA 



(5) 



where d is the length of the wall {x{s),y{s)), A the area 
of the positive domain (at one side of the wall, and 



\/8eoc 



(6) 



is the energy of the wall per unit of length. Consequently, 
the wall behaves as an elastic line with a linear tension 
a and pushed by a field H. 

Our approximations are exact for infinitely narrow in- 
terfaces, w — >■ 0, since we are reducing the field in the 
whole plane to a single curve defining the center of the 
interface. For thin interfaces the approximation is good 
enough, provided the width of the wall remains approx- 
imately constant all along the curve and that the local 
curvature radius of the line is smaller than the domain 
wall width w. In brief, these conditions assure that the 
state of (j) with a domain wall can be well described exclu- 
sively by the transverse degrees of freedom of an elastic 
interface. 



B. Interface pinning in a holed medium 

We will now consider a domain wall in a two dimen- 
sional medium with holes or multiply connected space. 
For the scalar field (p this amounts to solving Eq. ([1} in 
a domain £7 — A , which includes all the two dimensional 
space ft, except the possibly non-compact region A occu- 
pied by the holes and outer space. In order to model the 
absence of material in A , we set free (Neumann) bound- 
ary conditions for the order parameter, dn(f>\dA = at 
the hole borders and sample boundaries 9A. 



Figure 1: (color online) Interface in a two dimensional 
medium with free holes and sample boundaries. The gain 
of interface energy by optimally intersecting the holes and 
sample boundaries produces domain wall pinning. 



Figure 2: (color online) The energy of an elastic wall (red on- 
line) is minimized by an arc of radius r — a/{2H). An equilib- 
rium state is reached when the arc intersects the boundaries 
orthogonally. 



Within the interface approximation described in the 
previous subsection, interface pinning arises from the 
gain of line energy (reduction of the total length of the 
interface) that is possible by optimally intersecting the 
holes and sample boundaries (see Fig. [T]). From Eq. ^ 
the energy of the pinned domain wall then reads 



2HA 



(7) 



where di — /j*'^^ ds\/ -I- is the length of the inter- 
face segment connecting the holes i and (with i = 
and i = N designating the sample boundaries) and the 
area A now excludes regions belonging to A. 

The free (Neumann) boundary conditions for the order 
parameter at the hole and sample boundaries translate in 
the interface description in the orthogonality condition 



0, Vi 



(8) 



where is the tangent vector of the boundary and v.; = 
(i(si), 2/(si)) the tangent vector of the interface both at 
the intersection point {x{si),y{si)). 

Metastable states of the interface are therefore local 
minima of the energy ^ with segments satisfying the 
orthogonality constraint ([8]) at its ends. In the following 
we discuss the geometry of these optimal segments, which 
are the building blocks of our method. 



C. Equilibrium state of a wall 

Our next step is to calculate the equilibrium profile of 
an interface segment and its stability. The energy of an 
interface segment is both a function of its shape and the 
location of its ends or contact points. In order to find 
the possible metastable states of the segment we need to 



minimize the energy given by Eq. ([5]) with the constraint 

(ED. 

As shown in Appendix [B] the solution of the cor- 
responding Euler-Lagrange equation, regardless of any 
boundary condition, is a circular arc of radius 



2H 



(9) 



Consider now a wall confined between two irregular 
boundaries, as plotted in Fig. [21 Let / be the distance 
between the contact points, and 9i + 90° {82 + 90°) the 
angle formed by the upper (lower) boundary and the line 
connecting the two contact points. The elastic wall min- 
imizes its energy adopting the shape of an arc of radius 
r and it must be orthogonal to the boundaries at the 
contact points. As illustrated in Fig. [2j this implies that 
9, ^e-,^e and 



sin6' 



1/2 



HI 
a 



(10) 



where we have used the expression for the radius of the 
wall, Eq. 



D. Local depinning fields for an anchored wall 

We can now proceed and calculate local depinning 
fields for narrow domain walls bounded between two bor- 
ders, which are central for studying the ratchet effect. 
Given a metastable state of the anchored domain wall 
the local depinning field is defined as the maximum field 
it can support by deforming continuously as we increase 
the field. Above this local depinning field the domain 
wall escapes the local environment and slides until it is 
trapped again in a new metastable state with a larger de- 
pinning field, if it exists. Otherwise it continues sliding. 

As an illustration, consider the particular case where 
the bottom border is the x-axis and the top border 
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is given by and arbitrary smooth function f{x) > 0. 
The wall, as we have seen above, is an arc of radius 
r = <t/{2H). Its center must lie in the x-axis, say at xq, 
since the wall is perpendicular to the a:-axis at the lower 
contact point. The upper contact point, be- 
longs to the arc, hence {xi — xq)^ + f{xiY — r^, and the 
orthogonality condition implies: 



/(^i) 



Xi - Xq 

The upper contact point is then given by the condition 



(11) 



Note that only the solutions with f'{xi) > must be 
taken ii H > 0. If such a solution xi exists for a given r 
(i.e., for a given field H), the lower contact point is given 
by X2 = xi if f'{xi) = and X2 ^ r + xi - f{xi)/f{xi) 
otherwise. 

For _ff = (r — >■ oo) the only possible solutions are 
points xi such that f'{xi) = 0. These solutions are 
straight vertical segments joining the two borders at xi. 
For concreteness let us assume that xi = for H = and 
that /"(O) > 0, so the initial state is metastable. In such 
a case if we quasistatically increase H (decrease r) from 
xi = we can generate a continuum set of solutions xi (r) 
parametrized by the field. At some field He = cr/^rc it is 
possible however to have a discontinuity in xi{r) due to 
the absence of solutions beyond He ■ We can then define 
the critical radius of the initial metastable state as 



XI \ f'[Xi) 



(12) 



/'(^i)>0 



where the condition f'{xi) > assures that is positive, 
so we obtain the forward depinning field. The depinning 
field is therefore He = ajlrc-, the upper contact point 



of the critical arc is 

„low 



= xi(r'c) and the lower contact 
= r + Xc- f{xc)/f'{xc)- 
As an illustration consider the geometry displayed in 
Fig- El /(a;) = 2 — cos(x). As a fimction of the field 
H stable arcs have a radius r = a/2H and the upper 
contact point must satisfy 



2 - cos(a;) n . , ,„ 

sin(a;) 



(13) 



The solutions of this equation for different fields are 
shown graphically in Fig. [3] (upper plot), and the corre- 
sponding arcs are shown in Fig. |3] (lower plot). The criti- 
cal state (also shown in Fig.|3](lower plot)) has a;"P ~ 0.85 
and corresponds to Vc ~ 2.23 and ^2°^ « 1.30. The ini- 
tially flat interface for H = will shift forward quasistat- 
ically upon increasing the field, following the xi{r) curve. 
Above He — a/{2rc) the interface will move at a finite 
speed. 

Let us now analyze the depinning from a rounded tip 
of curvature radius VF, as shown in Fig. 21 We assume. 
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Figure 3: Construction of the metastable elastic lines bounded 
between and f{x) — 2 — cos{x), as a function of the ex- 
ternal field or arc radius. We depict the boundaries and the 
metastable walls for different fields H in the lower plot. In the 
upper plot we represent the function [f{x) //' (a;)] -^1 -I- f'{x)'^. 
The intersection of the arc radii r = a/{2H) with this function 
gives the upper contact point (a:i,/(a;i)) of the metastable 
wall with the top border f{x). The depinning field is He, 
above which there are no metastable states. 



for concreteness, the form 



fix) = {lo/2 + W)~Vw 
f{x) = f{xo) -|-tana|a; 



a;2 ,|a;| < xq (14) 
xo\ , |a;| > xo (15) 



where xq — VFsin(a) and f{xo) ~ (^o/2 + W) — Videos a. 
The first equation describes a rounded circular point, and 
the second a line with the asymptotic slope angle a. In 
this case xi increases monotonically from zero and no 
more solutions of Eq. (fTTj) exist for x > xq. We thus 
have Xc = VFsina and Vc = /(xq)/ sin a. The critical 
field is therefore 



Hr 



asm a 



lo + 2W{1 ~ cos a) 
For a sharp W ^ tip we have 



Hr 



asm a 



lo 



(16) 



(17) 



implying a very strong pinning in the limit of strong con- 
striction ^0 — 0. Interestingly, in this limit the depinning 
field would be ultimately controlled by the rounding W 
in more realistic rounded tips. 
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Figure 4: A tip rounded at the scale W . 




Figure 5: (color online) A two dimensional ratchet geometry 
made by triangles. 
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Figure 6: (color online) Broken symmetry in a 2D kinked 
wall, (a) Kink moving up. (b) After an x-axis reflection, an 
antikink moving down is obtained, (c) After a second y-axis 
reflection, a kink moving down on a different array (inverted 
triangle array) is obtained. From (d) to (f) same reflections 
are shown for an antikink. 



III. BUILDING A 2D RATCHET 



Prom this general theory of interface pinning, we will 
show how to build a 2D ratchet for extended domain 
walls with both direct and inverted rectification effects, 
as a function of the applied field, by choosing the appro- 
priate geometry for an array of asymmetric holes. We 
take as a starting point the geometry depicted in Fig. [SI 
where triangular defects are distributed in a rectangular 
array, which is similar to the hole arrangement analyzed 
in Refj^. This array of triangles presents a symmetry of 
reflection along the X axis but broken reflection symme- 
try along the Y axis, which is the basic condition for the 
observation of ratchet effects: the equivalence between 
forward and backward domain wall propagation is bro- 
ken in the array allowing for domain wall rectification 
effects. 

This kind of ratchet effect has been mostly studied in 
ID magnetic nanowiresi^— . However, in a 2D array 
such as shown in Fig. [S] the ID character of the elastic 
domain walls opens the possibility of extra propagation 



modes. In particular, when a wall is pinned between two 
lines of defects, it develops kinks and antikinks, as shown 
in Fig. ini Depending on their shape and the sign of the 
field, kinks and antikinks can move upward or downward, 
turning in a net wall motion to the left or to the right. 
This kink motion is also asymmetric, reflecting the Y- 
axis asymmetry of the pinning potential by the array of 
triangles, so that it opens the possibility of a rectified mo- 
tion of the kinked wall. In particular, let us analyze in 
detail how the symmetry properties of the array influence 
kink propagation: for example, as shown in Fig. [51 a kink 
moving upward (Fig. |6l^a)) is equivalent to an antikink 
moving downward (Fig. EJb)) upon reflection along the 
X-axis, which is an allowed symmetry operation of the ar- 
ray of triangles. In fact, both movements (kink upward 
and antikink downward) result in a net backward mo- 
tion of the extended domain wall. The critical field for 
this propagation process will be labelled H\j from now 
on. However, a kink moving downward (Fig. I^l^c)) is the 
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result of a reflection upon the Y axis of the antikink mov- 
ing downward. This is a broken symmetry in the array 
(note the inverted triangles), implying that both situa- 
tions are not equivalent. Actually, both the downward 
motion of a kink (Fig. [6je)) and the upward motion of 
an antikink (Fig. |6l^d)), which are equivalent upon reflec- 
tion along the X-axis, result in a net forward motion of 
the extended wall. The critical field for this propagation 
process will be labelled Hjy in the following. In short, 
H\j and Hd could not be the same due to the broken Y 
axis symmetry in the array. 

Thus, to understand domain wall propagation in the 
rectangular array of triangles two facts must be consid- 
ered: first, the broken Y-axis symmetry breaks the equiv- 
alence between forward and backward domain wall prop- 
agation; second, the extended nature of domain walls in 
the 2D array of holes allows for extra propagation modes 
not possible in ID geometries, such as those correspond- 
ing to nanowires: flat wall propagation and kinked wall 
propagation. If a wall is pushed to the right by the ap- 
plied field two different propagation modes can be acti- 
vated: either forward flat wall propagation (at a critical 
field Hp) or kinked wall propagation (by kinks moving 
downward and/or antikinks moving upward at Hd). On 
the contrary, if a wall is pushed to the left by the ap- 
plied field the possible propagation modes will be either 
backward flat wall propagation (at a critical field Hb) or 
kinked wall propagation (by kinks moving upward and / or 
antikinks moving downward at Hu). The global behavior 
of walls upon propagation across the array of asymmetric 
defects will depend on the relationships between the four 
relevant critical fields Hp, Hb, Hu and i?D- Actually, 
for some geometries, it can happen that the rectification 
of kink motion is opposite to the rectification of a ver- 
tical non-kinked wall described previously. This crossed 
ratchet effect offers promising technological applications 
since it allows a non trivial control of the two dimensional 
wall. 

First, we will analyze the propagation of a flat wall 
crossing a line of triangular defects (forward-backward 
ratchet), which is equivalent to the propagation of do- 
main walls in nanowires with asymmetric geometryi^Ti^^. 
Then, we will study the upward-downward propagation 
of a kink in a wall pinned in between two adjacent defect 
lines (upward-downward ratchet). Finally, we will discuss 
the geometrical parameters needed to design 2D arrays 
of asymmetric holes with crossed ratchet behavior (op- 
posite sign for forward-backward and upward-downward 
ratchets). 



A. Flat walls: Forward-Backward propagation 

An infinitely narrow domain wall moving from the left 
to the right (forward) across a vertical line of triangles 
(see Fig. [S]) will be pinned at the base of the triangles, 
where the distance between the ends of the wall is min- 
imum. The depinning field can be derived from the one 
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Figure 7: (color online) Scaled Hb/w vs. {Iq + w)/(Iq + b) 
for different values of c, lo and b. Solid line indicates the 
Hb/w ~ {lo + w)'^ dependence. Inset shows backward de- 
pinning field for a 0^ domain wall, Hb, as a function of the 
vertical gap between triangles, lo, for different values of the 
elastic constant c or domain wall width w. 



obtained in subsection lll DI for a geometry defined tip and 
a flat boundary, as in Fig. |4l Now we have two symmetric 
tips but this situation is equivalent to the previous one if 
we add the mirror image of the tip, the wall adopting the 
same shape as the one depicted in Fig. [4] (plus its mirror 
image). Therefore, the depinning field is identical to the 
one given by Eq. (fT7)) : 



Hp = 



(18) 



where 6 is the angle between the sides of the triangles 
and the horizontal, and lo is the vertical distance between 
triangles (see Fig. [5]). Notice also that the critical field 
(UHl) is, according to Eq. ([TU)) . the one for which the wall 
accommodates to the boundary of both triangles. 

A wall moving to the left (backward), like the one de- 
picted in Fig. O will be also pinned between the same 
vertices of the triangles, but now it has to grow along 
the vertical bases, i.e., the angle of the boundaries in 
Eq. (ini) is 90°. Therefore: 



Hb = ^ 
in 



(19) 



We find Hf/Hb = sin^ < 1, i.e., it is easier for the 
wall to move forward than backward, as expected. The 
triangles can therefore rectify the motion of the wall. 
Applying an alternating field of peak intensity H, with 
Hb < H < Hp, the wall will have a net forward motion, 
so that a direct ratchet effect is obtained. 

Finite width efTects. We compare our simple pre- 
vious geometrical estimates for Hp and Hb with simu- 
lations for a more realistic cj)'^ domain wall with a finite 
width, as it was done in Ref. S]. In the inset of Fig. [7] we 
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Figure 8: (color online) Forward/backward asymmetry for 
the flat 0'* wall with c = 50, Hf/Hb vs. sin6', calculated for 
different lo/{lo + b). Inset shows the forward depinning field, 
Hp for a (j)'^ domain wall, as a function of the isosceles angle, 
9, for different values of the vertical gap between triangles, Iq. 




Figure 9: (color online) Upward motion of a kink in an elastic 
wall (red online). 



plot the backward depinning field, Hb, of a 4>^ domain 
wall as a function of Iq for different elastic constants c 
and fixed eo = 1. This is equivalent to change the do- 
main wall width w — ■y/2c/eo and the domain wall energy 
a — v^Sceo/S that scales as a ^ w. Hb is found to in- 
crease as a function of c, mainly due to the increase in do- 
main wall energy a and to decrease as a function of ^o- All 
the data of Hb obtained from the different simulations 
can be scaled to a single curve, if we plot Hb/w ^ Hb/cf 
as a function of {lo+w)/{lo+b) (see main panel of Fig. [7]). 
For small values of la, Hb/w ^ (Iq + w)^^ as predicted 
by Eq. ([T9| except for a correction to Iq which is of the 
order of the domain wall width w. This correction can 
be qualitatively understood by noting that the (j)'^ wall 
just below Hb extends from the tip of the triangles up 
to a distance of order w into the base of each of them, 
so that the center of the wall describes an arc covering a 
vertical distance Iq + w. It is interesting to mention that 
this correction is the same as predicted by Eq. p6)) for 
a rounded tip of curvature radius equal to domain wall 
width {W — w) and a = 90° (as corresponds for back- 
ward depinning). That is, finite domain wall width and 
tip rounding of defect geometry have equivalent effects 
on depinning fields, softening the magnetic behavior in 
comparison to analytical calculation for sharp tips and 
narrow domain walls. At large values of Iq, Hb/w de- 
viates from the behavior Hb/w ^ {Iq + w)^^ decreasing 
steeply precisely when (Zq +w)/{lo + b) ^ 1. The reason 
is that for the simulations with varying Iq we fix the pe- 
riodicity of the lattice Iq + b. Therefore, for large values 
of Iq the base of triangles becomes small, and eventu- 
ally of order w, strongly reducing the geometric pinning 
mechanism when w ^ b. 

In the inset of Fig. [5]we show Hp vs. 9 for the 4>^ wall. 



for different values of Iq and constant c = 50. The main 
panel of Fig. [51 shows the forward/backward asymmetry 
for the flat wall calculated as Hp / Hb in comparison with 
the analytical prediction of Eq. (llSp for a narrow domain 
wall Hp/Hb = sin 6*. The simulated values follow nicely 
the sin 9 line except for small deviations at small and 
large angles 9 ^ 90°. These can be in part attributed to 
the discreteness of the lattice, which does not allow to 
produce smooth slopes at the scale w when the angle is 
too close to 6* = and 6* = 90°. 

In short, our simulations with the 0^ model are consis- 
tent with the geometric estimates for a narrow wall and 
show how to correct the depinning fields for single arcs 
with a finite width, which can be relevant for experimen- 
tal situationsSi,. The depinning of single arcs are, on the 
other hand, the main building blocks for calculating all 
the anisotropic depinning fields and, in particular, the 
crossed-ratchet effect. Thus, the softening of the criti- 
cal fields observed due to finite width corrections and / or 
the effect of rounded tips could also be applied in a simi- 
lar way to the geometric estimates of the propagation of 
kinked walls. 



B. Kinked wall: Upward-Downward propagation 

Upward propagation. The depinning field for kinks 
can also be calculated using the basic geometry analyzed 
in Sec. Ill CI Fig. [9] shows the evolution of a wall forming 
a kink as it is pushed upwards when the field increases 
from H — 0. The initial disposition of the wall is labelled 
as 0. The critical field Hij = (T/(2rmin) for a complete 
depinning of the kink is given by the arc with minimal 
radius r,„i,i. For a given geometry, one has to carefully 
trace the trajectory of the wall, as depicted in Fig. [HI and 
compute the minimal radius in each step. 
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To step from to 1, i.e. to depin the transverse hor- 
izontal segment of the kink, is similar to the situation 
depicted in Fig. |4] with a = 90° — 0. Therefore, the corre- 
sponding critical field is ctcos0/(2/i). However, the base 
h of the triangle can be too short for the domain wall to 
reach the symmetrical position described in Fig. |4l The 
wall is maximally tilted at an angle /3 (see Fig. [5|). In 
this situation, the corresponding angles d\ = 62 in Fig. [5] 
are equal to /? (the angle formed with the base of the 
rightmost triangle in Fig. |9]) and the distance between 
the two contact points \s I = h/ cos/3. The wall reaches 
this orientation \i H > asmP/l = crsin(2/3)/(2/i). Con- 
sequently, the critical field to move from to 1 is the 
minimum of these two fields, namely. 



a 
2h 

|-minjcos(0),^^ 



min {cos(0), sin(2/3)} 
bh 



{b/2y 



(20) 



The next critical arc is 3, a wall orthogonal to the 
contact sides of the triangles across a diagonal of the 
rectangular cell of triangles. The center of this arc is the 
point C3, the intersection between the two prolongations 
of the triangles sides. The radius of arc 3 is + b/2 + 
/i/tan0, and the corresponding critical field: 



1 



2h V ^o/Zi + tan /3 + 1/ tan 9 



(21) 



Finally, one should also consider the diagonal arc 5 to 
complete the upward motion of the kink. However, the 
radius of this arc is bigger than ht + h — 6/(2 tan0) -I- h, 
resulting in lower critical fields for the geometries con- 
sidered in this paper. The final result for upward motion 
is: 



(22) 



From the above equations, three geometrical param- 
eters of the rectangular array of triangles are found to 
control de interplay between the different depinning pro- 
cesses of the kinked wall and, thus, the relevant critical 
fields: the angle that defines triangle shape, the angle 
(3 that characterizes the shape of the horizontal inter- 
triangle region (/3 is given by tan/3 = b/2h), and the 
ratio h/lo between horizontal and vertical triangle dis- 
tance. This last parameter, h/lo, is only important in 
the depinning of the diagonal arc 3. Figure [TU] shows the 
calculated Htj, normalized by the scale factor a/2h, as a 
function of /3 for = 45° and different values of the ratio 
h/lo — 10,1,0.1. For large /3, Hu is given by , so 
that it is the same in the three panels of Fig. (TU] Below 
P ~ 22°, there is a crossover to Htj — H^ indicated by 
the upturn in i?u(/3) as /3 decreases. It occurs at different 
angular positions depending on h/lo- Pc = 2.6°, 13° and 
21° for h/lo = 0.1, 1, 10, respectively. That is, for small /3 
and large h/lo (very anisotropic rectangular array), crit- 
ical upward depinning occurs at the diagonal arc 3 in 
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Figure 10: Critical fields, and H\j normalized by cr/2/i vs. 
P calculated for 6 = 45° for (a) H/Iq = 10 (b) h/lo = 1 and 
(c) h/lo ~ 0.1. Note that for large /3, Hu becomes equal to 
Hu and, therefore, a single line appears in the figures. 



Fig. ini whereas in the rest of the parameter space the 
relevant process corresponds to depinning of the trans- 
verse domain wall segment between adjacent triangles in 
the same row (arc 1 in Fig. [9]). 

Downward propagation. The first step — ^ 1 of 
the downward motion (Fig. [TT|) is identical to the same 
step in the upward motion, hence Hf = H^ ■ From 1 
to 3, the arc 2 has the minimal radius, which equals the 
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Figure 11: (color online) Downward motion of a kink in an 
elastic wall (red online). 



distance between its center A and the vertex Vi, {h + 
ht)sm9. However, point B can lie below or above the 
side of the triangle. In the first case, which occurs if 
> 45°, the critical field is given by the radius of arc 3, 
(h+ht) / {2 cos 9). The second case occurs when < 45° — 
/3 and then the minimum radius, h/{2cos{(3 + 6) cos (3) 
occurs when point B is at the vertex of the triangle. The 
field to reach arc 3 can be written as: 



^ 2h 



( 2cos(/3 + 6')cos/3 
1 

[tan/3/tan6'+ 1] sin 9 
2 cos 6* 



if 6* < 45° - /3 



if 45° - /3 < 6* < 45° 



if 61 > 45° 



I tan /3/ tan 6* + 1 
Finally, the radius of the diagonal arc 4 is 




2lob + b^\ 



yielding 



a 

2^ 



(23) 



(24) 



different values of the ratio h/lo = 10,1,0.1. For large 
/3, Hy) is given by Hf, i.e. depinning of the transverse 
horizontal wall segment, but as /3 decreases and Hf 
become more relevant. In particular, for large h/la, 
dominates the behavior in a wide (3 range. The result 
is that, in the low /3 range, Hd is much larger than Hi], 
but above a certain threshold Pq both fields become equal 
{Hb = Hij) (for example, for h/lo = 1, = 45°). Thus, 
for /3 < /?o , upward kink propagation is easier than down- 
ward propagation so that when an alternating field of 
peak intensity H, with Hu < H < Hiy is applied to 
the kinked wall, it will have a net backward motion, i.e. 
opposite to the behavior observed in the previous subsec- 
tion on flat wall propagation walls. On the other hand, 
for /3 > /3q, Hd = Hu, so that kink propagation is not 
rectified by the array of triangles. 

In short, for 9 = 45°, as is the case in FigfTUl whenever 
kinked wall propagation is asymmetric, it results in an in- 
verted ratchet effect. This is actually the case for most of 
the parameter space (/3, 9, H/Iq). For example, Fig.[T2la) 
shows the phase diagram in the (/3, 9) plane for H/Iq — 1, 
in which only these two regimes for domain wall propaga- 
tion are found: inverted ratchet {Hd > H\j) in the low /3 
region and symmetric kink propagation in the right bot- 
tom corner of the diagram. This is a direct consequence 
of the maximum condition imposed in eq. (j26p . as long 
as H^ = Hi' . However, at large h/lo, the role of arc 3 
in Fig. [9] in critical upward depinning becomes more im- 
portant and H^ is given by H^ in a wider (/?, 9) region. 
In this case, Hu can take any value in comparison with 
Hd , so that a direct ratchet effect for kinked wall motion 
becomes possible. An example of this situation can be 
seen in Fig. HHb) for 9 = 78° and H/Iq = 10. Thus, the 
phase diagram for this very anisotropic array of triangles 
with h/lo = 10 (Fig. [12] (c)) becomes more complex: in- 
verse ratchet effect (i/o > ^^u) is found in a large {(3,9) 
region in the left part of the diagram, kink motion is 
symmetric in the right part of the diagram (iJn = ^u) 
and, finally, direct ratchet (i.e. Hd < Hu) is found in 
two small regions close to the upper part of the diagram, 
above 9 = 65°. This direct ratchet regions shrink as h/lo 
decreases and disappear for h/lo < 2 due to the soft- 
ening of the depinning processes of the diagonal arcs in 
comparison with depinning of the horizontal transverse 
segments. 



Crossed ratchets 



and 



'2h [y'ilo/hy + + lo/h tan /? + tan^ 



=max{ffU,iJ3°,fff} 



(25) 



(26) 



Figure [TU] shows the calculated H-q , normalized by the 
scale factor a/2h, as a function of /3 ior 9 — 45° and 



From the previous analysis, we have found a fundamen- 
tal difference between flat and kinked wall propagation 
modes: H-p/H^ is always smaller than unity, implying 
that flat wall propagation under an ac field will result in 
direct ratchet effects; on the contrary, Hu/Hu, can take 
any value so that kinked wall motion can result either 
in direct and inverse rectification effects. Thus, the first 
condition to design an asymmetric array of defects that 
displays crossed ratchet behavior is to choose a point in 
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Figure 12: (color online) (a) Phase diagram in 6') plane 
for of the different regimes for kinked wall propagation for 
h/lo — 1; (b) Critical fields, Ho and H\j normalized by a/2h 
vs. P calculated for 6 = 78° for h/lo = 10; (c) Same as in 
(a) for h/lo = 10. Dotted lines in (a) and (c) indicate the 
condition Ho ~ Hp below which flat wall propagation modes 
compete with kinked wall propagation. 



the phase diagram of Fig.lT^in which Ht]/Hn < 1. Then, 
in order to observe clear crossed ratchet effects that can 
be useful for device applications, the interplay between 
the four relevant critical fields Hy , Hb , -ffn and H\] must 
be taken into account. 

In Fig. [13] we plot the four critical fields as a func- 




Figure 13: (color online) Critical fields, for 6 = 30°, and 
P = 30° as a function of h/ lo- 



tion of h/lo between triangles, normalized by for 
(3 — 30° and 9 = 30°, which is similar to the geometry 
used in the experiments of Ref. |3i, where crossed ratchet 
effects were observed both experimentally and theoreti- 
cally. The crossed ratchet effect is apparent from the fig- 
ure: Hp is smaller than Hq, but Hjj (the upward motion 
of the kink drives the wall backwards) is larger than Hd. 
Domain wall propagation in the array is determined by 
the relationships between the four critical fields: for each 
particular array geometry (i.e. each H/Iq value) when 
a domain wall is pushed in the forward direction it will 
depin as soon as the applied field reaches the lowest of 
Hp and Hd; but, when a domain wall is pushed in the 
backward direction it will depin as soon as the applied 
field reaches the lowest of Hb and H\j. For large h/lo 
(i.e. very rectangular array cell), both Hd and Hu take 
much lower values than Hb and Hp implying that two 
well separated field ranges can be defined: low field do- 
main wall propagation dominated by kink motion (i.e. 
easier backward wall motion) and high field domain wall 
propagation dominated by flat wall motion (i.e. easier 
forward wall motion). As h/lo is reduced below w 1, Hd 
becomes larger than Hp, the interplay between flat and 
kinked wall propagation modes becomes more complex, 
and the different rectification effects cannot be clearly 
separated. Finally, for very close triangle lines {h/lo be- 
low 0.25), Hu becomes larger than Hb and domain wall 
motion in the array is dominated by flat wall propagation 
modes. 

The condition Hd = Hp is plotted as a dotted line 
in Figs. [12] (a) and (c), so that the region for well sep- 
arated kinked and flat propagation modes, i.e. clear 
crossed ratchet observation, lies above this line in the 
{13,6) plane. It can be seen that as h/lo increases the 
available parameter region for crossed ratchet becomes 
wider due to the different scaling of the critical fields: 
Hp and Hb scale as I/Iq, whereas Hd and Hu scale as 
l/h. Thus, the design of arrays in the large h range 
{h ^ ^o) appears as an important condition for a clear 
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observation of crossed ratchet effects that can be of use 
in device apphcations. 



IV. CONCLUSIONS 

In summary, the propagation of an elastic domain wall 
in a two dimensional medium has been analyzed in an 
arbitrary geometry defined by holes and sample bound- 
aries. The local depinning fields for an anchored wall 
have been calculated as a function of boundary shape in 
terms of the minimal arc radius that satisfies the rele- 
vant orthogonality conditions. Then, these results have 
been applied to the design of 2D arrays of asymmet- 
ric holes with broken Y-axis reflection symmetry that 
can display crossed ratchet effects (i.e. direct ratchet 
for forward/backward flat wall propagation and inverted 
ratchet for upward/downward kink propagation). 

For a rectangular array of triangles, fiat wall propaga- 
tion is found to be always asymmetric resulting in a di- 
rect ratchet effect controlled by triangle shape (angle 9) 
and intertriangle vertical distance (Iq). Corrections due 
to finite domain wall width and/or rounded triangle tips, 
that could be relevant in actual patterned arrays of holes, 
produce a global softening of the critical fields but do 
not alter significantly forward/backward asymmetry. On 
the other hand, upward/downward kink propagation can 
display any asymmetry and depends not only on triangle 
shape but also on the shape of the horizonal intertrian- 
gle region (angle /?) and on the array vertical/horizontal 
anisotropy {h/lo). The array geometry needed for the ob- 
servation of crossed ratchet effects has been determined 
considering the different wall propagation modes relevant 
in the different points of the (/3, 9) plane. Anisotropic 
arrays with large h/l^ are found to be optimum for the 
observation of clear crossed ratchet effects. 
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Appendix A: From field equations to elastic 
interfaces 

In this Appendix we construct an approximate station- 
ary solution of Eq. ^ around a given curve in the plane. 



defined as {x{s), y{s)), with s a real number taking values 
in some interval. 

The interface of the desired solution is centered around 
the line (x(s),y(s)), i.e., (f)(x{s),y(s)) — 0, and the field 
approaches to the stable values ±1 as we move away from 
the line. We then construct the solution using the signed 
distance function, g{x, y), whose absolute value is the dis- 
tance of a point {x,y) to the line {x{s),y{s)). Obviously, 
g{x{s),y{s)) = for all s. One less obvious property 
is that the gradient of the distance function is unitary 
all over the plane. In other words, the distance function 
obeys the eikonal equation: 



[dxg{x,y)f + [dyg{x,y)f = 1 



(Al) 



Now we choose the following form for the field 
(j){x,y) — f{g{x,y)). Introducing this ansatz in the sta- 
tionary ip"^ equation and making use of Eq. (|Aip . we get: 

cf'ig) + cf'{g)y^g + e,[f{g) - fig))'] = (A2) 

Our first approximation consist of neglecting 
f'{g{x,y))V'^g{x,y) in the above equation. The 
Laplacian of the distance function is inversely propor- 
tional to the curvature radius of the line {x{s),y{s)). 
Therefore, our approximation is valid for interfaces with 
a curvature radius much larger than its width. Then, 
Eq. (jA2l) reduces to: 



c/" (g(x, y)) + eo[f{g{x, yj) - f{g{x, y))^] - (A3) 

and the general solution reads f(z) = tanh[(z — zo)/w] 
with 



w 



The field (j) is then given by: 



(x, y) — tanh 



.9(2;, y) 



(A4) 



(A5) 



where we have absorbed the constant zq in the function g 
to center the wall along the line (x(s), y{s)) where g van- 
ishes. Eq. (jASP is well known as the field corresponding 
to a single wall in the one dimensional (j)'^ model. 

To calculate the energy of the solution given by 
Eq. (|A5p . it is convenient to use as coordinates the dis- 
tance z to the center of the interface and s, the pa- 
rameter defining this center. These new coordinates 
(s, z) are related with the cartesian coordinates (x, y) as 
{x{s,z),y{s,z)), obeying: 



g{x{s,z),y{s,z)) ^ z ys,z 



(A6) 



The z — contour line is our initial curve {x{s),y{s)). 
The Jacobian of this change of coordinates can be cal- 
culated differentiating Eq. (|A6I) with respect to s and z, 
respectively, yielding 



dxdy = \/ x'^ + y'^ds dz 



(A7) 
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where the dot denotes differentiation with respect to s. 
This is in fact the product of the elementary length of 
the contour line y'x^+^ds times dz. 

We can now calculate the energy of the wall inserting 
the solution (lASP in Eq. ([2]). With the change of variable 
(x, y) — ?> (s, z), the energy reduces to: 



E= I dzL{z) U(f{z)) + Hf(z) + -f{zf 



where f{z) ~ tanh(z/w) and 



L{z) — I ds\/ i{s, z)2 + y{s, z)^ 



(A8) 



(A9) 



is the total length of the contour line g{x,y) = z in the 
restricted geometry of the problem. In particular, L(0) 
is the length of our starting curve {x{s),y{s)) defining 
the center of the interface. If the interface is narrow, we 
can approximate L{z) ~ ^(0) for those z where U{f{z)) 
is significantly different from zero, i.e., around the center 
of the interface. Finally, the energy due to the external 
field H can be estimated replacing /(z) by a step function 
26'(z) — 1 in the term Hf{z). With these assumptions, 
the energy becomes: 



E = crL(O) - 2HA 



(AlO) 



where A is the area at one side of the center of the inter- 
face (x(s),y(s)), and 



V8eoc 



(All) 



is the energy of the interface per unit of length. These 
are the expressions yielding ([5]) in the main text. 



Appendix B: Equilibrium shape 

Our next step is to calculate the equilibrium profile of 
an interface segment. The energy of an interface segment 
is both a function of its shape and the location of its ends 
or contact points. In order to find the possible metastable 
states of the segment we need to minimize this energy 
with the constraint ([5]). 

We describe the wall as a line given by y{x), anchored 
to the border at points {xi,yi) on the left and (x2,?/2) 
on the right (see Fig. [T4)) . The energy of the interface is 
given by 



E : 



a^/TT7W-2Hy{x) 



dx. 



(Bl) 



hence the Euler-Lagrange equation is 



d y' 2H 



dx ^/mw <y 



= 0, 



(B2) 




Figure 14: (color online) The elastic wall (red online) between 
two boundaries (vertical black curves) is parametrized as y{x) 
to solve the variational problem. 



We have to solve this equation imposing the orthogonal- 
ity condition at the contact points (which are otherwise 
free). One integration of (|B2|) gives 



X — Xo 



(B3) 



where Xq is a constant and r = a/{2H). From (jB3[) we 
get 



y'{x) = ±- 



X — Xq 



and a second integration yields 



(B4) 



y{x) ± - (a; - x^f 



(B5) 



which, written as (x— a:o)^-l-(?/— yo)^ = ''^ reveals itself as 
the arc of a circumference of radius r and center (a;o, yo)- 
This is in fact the Laplace law in two dimensions, relat- 
ing the pressure difference to the local curvature of an 
elastic interface at equilibrium. It is however important 
to notice that this equilibrium shape is independent of 
the boundary conditions. We can therefore impose these 
conditions by looking for an arc of radius r which inter- 
sect orthogonally with the two boundaries, as we do in 
Sec. Ill Cl using basic geometric arguments and in Sec. Ill Dl 
in an analytical manner. 
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